Effective temperature of an aging powder 
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The aging dynamics and the fluctuation-dissipation relation between the spontaneous diffusion 
induced by a random noise and the drift motion induced by a small stirring force are numerically 
investigated in a 3D schematic model of compacting powder: a gravity-driven lattice-gas with purely 
kinetic constraints. The compaction dynamics is characterized by a super-aging behaviour and, in 
analogy with glasses, exhibits a purely dynamical time-scale dependent effective temperature. A 
simple experiment to measure this quantity is suggested. 
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Introduction. - Slow relaxation phenomena are ubiqui- 
tous in nature. When one deals with glassy dynamics the 
challenge is the identification of the relevant degrees of 
freedom which makes a thermodynamic description still 
possible. In the attempt to provide a unifying framework 
to the behaviour of aging and non-relaxational systems a 
microscopic definition of effective temperature was intro- 
duced through a generalized fluctuation-dissipation rela- 
tion [0. This quantity turns out to coincide, at least in 
mean-field glasses, with the Edwards' compactivity M 5|, 
previously introduced in a granular matter context 
For a class of finite-dimensional and zero-gravity com- 
pacting systems, recent numerical results have come to 
support this correspondence However, in contrast to 
the case of strong ||, or moderately strong vibra- 
tion regime, the absence of a temperature-like quantity 
in slowly driven compacting systems under gravity is a 
rather problematic issue ]iT|-p^|. In this Rapid Com- 
munication we investigate the nature of the aging dy- 
namics and the effective temperature in a simple 3D 
lattice-gas model of powder |Q . We show that the non- 
equilibrium dynamics is characterized by a three-steps re- 
laxation mechanism and 'super-aging' in the mean-square 
displacement. We then find that during the compaction 
the response to a random perturbation is positive and 
observe a violation of the fluctuation-dissipation relation 
similar to glasses. Such features show that it is possible 
to describe the gravity-driven compaction dynamics in 
terms of a purely dynamical time-scale dependent effec- 
tive temperature. 



The model. - The model we consider was introduced 
in (13 as a simple generalization of the Kob-Andersen 
model jl5| . The system consists of N particles on a body 
centred cubic lattice where there can be at most one par- 
ticle per site. There is no cohesion energy among parti- 
cles and the Hamiltonian is simply 



N 



ticle i, and m its mass. At each time step a particle can 
move with probability p to a neighboring empty site if the 
particle has less than v nearest neighbors before and af- 
ter it has moved. Here p = min[l, x~ Ah ] where Ah = ±1 
is the vertical displacement in the attempted elementary 
move and x = exp(— mg/T) represents the 'vibration'. 
The kinetic rule is time-reversible and hence the detailed 
balance is satisfied. We therefore assume that statistical 
properties of the mechanical vibrations on the box can 
be described, after a suitable coarse-graining, as a ther- 
mal bath at temperature T. In the regime of quasi-static 
flow this assumption, which neglects the complicated ef- 
fects like friction and dissipation between grains, is a rea- 
sonable starting point ||. We set throughout the mass 
to = 1 and the threshold v — 5. Particles are confined 
in a box closed at both ends and with periodic bound- 
ary condition in the horizontal direction. We consider 
a system of height 16L and transverse surface L 2 with 
L = 20, and number of particles N = 16000 (correspond- 
ing to a global density of 0.25). As shown in |l4| ], the 
interplay of kinetic constraints and gravity are enough to 
reproduce the basic aspect of weakly vibrated powder like 
compaction and segregation phenomena, and vibration- 
dependent asymptotic packing density. In the following 
we consider the non-equilibrium features as they show up 
in the two-times correlation and response function. 

Aging dynamics. - Slow relaxation in weakly vibrated 
powder is closely related to the reduction of the free vol- 
ume available to the particle motion, hence age depen- 
dent properties are expected |l?}-[i~8|] . The aging dynam- 
ics in simple models of granular matter has been recently 
studied by several authors ||l9|Jl2] , ^3| , po[ . For our pur- 
poses it can be easily characterized in terms of the 'mean- 
square displacement' between two configurations at time 
t,„ and t > L„: 



Hq = mg ^ h i 



(1) 



1 N 

B{t,t w ) = -J2{Mt) - Ht w )] 



where g is the gravity constant, hi is the height of the par- 



(2) 



where the angular brackets denote the average over the 
random noise. The system is initially prepared in a ran- 
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dom loose packed state, which in this model corresponds 
to a packing density p t \ v ~ 0.707 [Q, and then the vibra- 
tion is turned on. The plot of B(t, t w ), see Fig. [|a, clearly 
shows the well known aging effect. The system does not 
reach any equilibrium state on the observation time-scale, 
but rather persists in a non-stationary regime: the par- 
ticle displacements become slower and slower as the age 
of the system increases. Interestingly, the Fig. [l]a shows 
a three steps relaxation mechanism: a short-time normal 
diffusion; an intermediate sub-diffusive regime which is 
tempting to associate to the 'cage rearrangement'; and 
finally a relatively faster but still sub-diffusive regime 
which can be figured out ctS Sb C3j ge-diffusion'. 
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FIG. 1. a) Mean-square displacement B(t,t w ) vs t — t w 
measured from a random loose packed initial state 
(prip — 0.707). The vibration amplitude is x = 0.2, and 
waiting times are t w = 2 k with k — 10, 17 (from the top to 
the bottom). The dashed line represent the short-time nor- 
mal diffusion regime, b) The same data plotted vs the scaling 
variable — t 1_M . The super-aging exponent is /j, ~ 1.41 
for x = 0.2. 



Inspired by recent analytic results on a related anoma- 
lous diffusion model pl| , we use a scaling function rep- 
resenting a super-aging behaviour: in the Fig. [j]b the 
mean-square displacement data are plotted vs the vari- 
able t]^^ — t 1 "^. We find that the value fx ~ 1.41 gives a 
quite good data collapse for vibration x = 0.2 (Fig. |l|b). 
The exponent \i appears to increase with the vibration, at 
least on the accessible time-scale of our simulations. We 
obtain for instance /i ~ 1.28 for x = 0.05 and ll ~ 1.48 for 



x = 0.4. According the theory presented in |2l| the ex- 
ponent is universal and its value is fi ~ 1.48. We believe 
that the discrepancy observed at small x is due to the 
fact that the system is still far from the asymptotic high 
packing density regime where the predictions of ref. ]2l| ] 
apply. In the super-aging regime, the characteristic re- 
laxation time of the system grows faster than its age, a 
rather unusual feature as compared to glassy polymers 
or random magnets |l8| . 

Effective temperature. - The possibility of a thermody- 
namic description of slow relaxing systems is apparently 
ruled out by the breakdown of time-translation invari- 
ance, i.e. the presence of aging phenomena. In granular 
materials further complications may arise from the pres- 
ence of spatial inhomogeneities induced by the boundary 
conditions and the gravity direction. In analogy with 
glassy systems we show, however, that the violation of 
the fluctuation-dissipation relation is not arbitrary. It is 
such that the degrees of freedom associated with the slow 
motion can be considered as equilibrated at an effective 
temperature (vibration) higher than the one imposed by 
the external bath (forcing). In order to see this we need 
to compute the dynamical response function. We apply 
to the system a small random stirring force at time t w : 



N 

i=i 



(3) 



where /; = ±1 independently for each particle. The lin- 
ear regime is probed for small enough values of the per- 
turbation e. The integrated response function conjugated 
to Eq. (0) is the 'staggered displacement' 



1 N 

K(t,t w ) = -Y,(hMt) - h^t 



(4) 



where the overline denotes the average over the random 
stirring force. At thermal equilibrium k and B are time- 
translation invariant and the Einstein relation holds, 



K,(t t w ^ 



2T 



B(t - t v 



(5) 



It has been suggested that a simple generalization of the 
previous relation in the aging regime provides, in a suit- 
able long-time limit, a reliable definition of effective tem- 
perature fl: 



(G) 



In a class of solvable mean- field models of glasses [ ^7[ , 
T e g(t,t w ) has the following properties. When t — t w ~ 
0(t w ), T cS {t,t w ) = T; while for t/t w > 0(1), T cS (t,t w ) 
is a constant or slowly waiting-time dependent quantity 
higher than the bath temperature T. In recent exper- 
iments, effective temperatures have been measured in 
glycerol [E2| and laponite |23|]. 
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FIG. 2. Measuring the non-equilibrium fluctuation dissi- 
pation relation in a compaction experiment. The system is 
prepared in a random loose packed state, p r i p ~ 0.707, and 
then shaken with a vibration x — 0.2. The system com- 
pacts for a waiting time t w = 2 14 after which a small stir- 
ring force is turned on. a) Comparison of the mean-square 
displacement B(t,t w ) and the conjugated response function 
\(t,t w ) — 2T n(t,t w )/e vs t — t w . b) Parametric plot of 
B(t,t m ) vs \(t,t w ). The slope of the dashed line is ~ 0.22. 
The solid line is the equilibrium fluctuation-dissipation the- 
orem which is correctly recovered in the zero-gravity limit, 
when particles diffuse in the whole box and the system reaches 
a low-density regime. 

In order to measure the effective temperature the sys- 
tem is first prepared in a random-loose packed state. 
Then it compacts for a time t w under a given vibration, 
x. At t w the random perturbation is applied, and the 
staggered displacement between configurations at times 
t w and t is measured. Following this method we find a 
positive monotonic response function. At small time sep- 
aration the integrated response and correlation are equal 
while at later time they depart from each other, with the 
response increasing slower than the corresponding corre- 
lation (see Fig. ||a). The nature of the violation of the 
fluctuation-dissipation relation is best represented in the 
parametric plot of n{t, t w ) vs B{t, t w ). In Fig. |^b is shown 
such a plot for a waiting time t w = 2 14 MCs and vibra- 
tion x = 0.2 In close analogy with glassy systems we ob- 
serve two 'quasi-equilibrium' regimes. At short time the 
usual equilibrium fluctuation-dissipation relation holds 
and T e g(t,t w ) — T. At later time a strong violation of 



the fluctuation-dissipation relation is observed. In this 
regime the numerical data can be fitted with an excel- 
lent approximation by a straight line. This means that 
the slope of the parametric plot does not depend on the 
observation time and therefore the effective temperature 
may only depend on t w . A closer inspection of Figs. |l| 
and ^| reveals that the crossover between the two quasi- 
equilibrium regimes takes place over a time scale corre- 
sponding to the 'cage-rearrangement' motion. The pic- 
ture outlined above holds in the whole range of weak 
vibrations 0.05 < x < 0.4 and waiting times t w < 10'° 
that we have explored. In particular, we find that the ef- 
fective temperature decreases with the external vibration 
(bath temperature) and slowly with the waiting time t w . 

Finally, to clarify some problems raised in Jll[ , and dis- 
cussed in detail in |L2| , we have also studied the response 
function to a uniform, either positive or negative, force 
field. In both cases we find a non monotonic response 
function with a negative component. In particular, when 
kinetic constraints are removed and the system is at equi- 
librium, the thermal bath temperature is not recovered 
with this method. This shows the importance of using 
stochastic perturbations to define the temperature from 
the fluctuation-dissipation relation. 

Conclusion. - The notion of effective temperature has 
been often invoked as the first step towards the formu- 
lation of a non-equilibrium statistical mechanics for sys- 
tems as different as vibrated powders ||f7| , turbulent flu- 
ids p4| , and structural glasses We have shown in 
this paper that a purely dynamical time-scale dependent 
effective temperature appears during the compaction dy- 
namics of 3D constrained lattice-gas models of aging pow- 
der. Mean-field glassy models suggest that such a fea- 
ture is robust with respect to non-relaxational perturba- 
tions and time-dependent driving forces [B5] . Prelim- 
inary results confirm this expectation in more realistic, 
Hertz contact mechanics based models pq| . 

It is worth noting that the perturbed Hamiltonian (||) 
describes a particle system with a random distribution 
of mass with average m and standard deviation propor- 
tional to e. Hence the measure of the effective tem- 
perature presented here could be experimentally carried 
out in a suitably prepared sample of glass beads, pro- 
vided that e is small enough to probe the linear response 
regime, and that mass-segregation effects are negligible. 
With the particle tracking experimental facilities (such 
as, for instance, PEPT [^J) this measure might not be 
out of reach. 
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